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We calculate the resonant x-ray scattering (RXS) spectra near the Ce Lm absorption edge 
in CeBe, on the basis of a microscopic model that the 4/ states of Ce are atomic while the 5d 
states form an energy band with a reasonable density of states. In the initial state, we employ an 
effective Hamiltonian of Shiina et al. in the antiferro-quadrupole (AFQ) ordering phase, while we 
construct the wave function consistent with the neutron scattering experiment in the magnetic 
ground state. In the intermediate state, we take full account of the intra-atomic Coulomb 
interaction. Without assuming any lattice distortion, we obtain sufficient RXS intensities on 
the AFQ superlattice spot. We obtain the spectral shape, the temperature and magnetic field 
dependences in good agreement with the experiment, thus demonstrating the mechanism that 
the intensity is brought about by the modulation of 5ci states through the anisotropic term of 
the 5d-4/ Coulomb interaction. In the magnetic ground state, a small pre-edge peak is found 
by the E2 process. On the magnetic superlattice spot, we get a finite but considerably small 
intensity. The magnetic form factor is briefly discussed. 
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§1. Introduction 

Resonant x-ray scattering (RXS) has recently at- 
tracted much interest as a useful tool to investigate the 
orbital order, which neutron scattering experiments are 
usually difficult to probe. The resonant enhancement for 
the prohibited Bragg reflection corresponding to the or- 
bital order has been observed in several transition-metal 
compounds by using synchrotron radiaiiimpSYith photon 
energy around the K absorption edgcBa'crBu' 

For such K-edge resonances, Ap states of transition 
metals are involved in the intermediate state in the elec- 
tric dipolar {Ei) process, and they have to be modu- 
lated in accordance with the orbital order for observing 
signals. At the early statge, for LaMnOa, such a modu- 
lation was considered to come from the anisotropic, term 
of the Ap-id intra-atomic Coulomb interaction,EHZl^ but 
subseauepl jStudies based on the band structure calcu- 
latiorJa-BliEfO have revealed that the modulation comes 
mainly from the lattice distortion via the oxygen poten- 
tial on the neighboring sites. Similar conclusions haizje 
been obtained for t2g-electron systems, such as YTiOatj^ 
and YVO3.II3' This is because Ap states are so extend- 
ing in space that they are very sensitive to the electronic 
structure at neighboring sites. 

Not only transition-metal compounds but also rare- 
earth-metal compounds show the orbital order (usually 
an ordering of quadrupole moments). Recently, RXS 
experiments were carried out around the Ce Lm ab- 
sorption edge in CeBg, in which resonant enhancements 
have beerviQUfid on the quadrupolar ordering superlat- 



tice spots.li3'llj) In particular, Nakao et al^ have found 
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a one-peak structure as a function of photon energy on 
the spot G — (555), which was assigned to the Ei pro- 
cess. Here G is the scattering vector in units of 27r/a with 
a being the lattice constant. They have also measured 
the temperature and the magnetic field dependences of 
the intensities. The purpose of this paper is to analyze 
their experimental result on the basis of a microscopic 
model and thereby to elucidate the mechanism for RXS 
in CeBg. Some of the results-jreported here were briefly 
presented in a recent letter.113' In this paper, making a 
slight revision on the bd density of states, we describe 
explicitly the model as well as the calculational proce- 
dure. We also add the calculation of the RXS spectra on 
the magnetic ground state. 

Each Ce atom is considered to be in the f^- 
configuration, ^-^5/2. The Fg quartet states have a lower 
energy than the F7 doublet under the cubic crystal field. 
The 4/ states are assumed to be atomic as a first approx- 
imation. With decreasing temperatures the antiferro- 
quadrupole (AFQ) order appears at Tq — 3.2 K with an 
ordering wave vector Q — (^5^), as shown in Fig. 1(a). 
This phase transition originates from the intersite inter- 
action between the atomic Fg states by lifting the degen- 
eracy of the Fg states. Ohkawa derived an effective pfp;- 



site interaction on theii 
Recently, Shiina et "''^ 



^RKKY interactionO'EJ' 
£3 extended his model by 



taking full account of the symmetry of the interaction 
as well as the order parameters. Thereby, they solved 
a lopscstanding ccptroversy between the neutron diffrac- 
tioncj) and NMRBS' in the context of the induced order 
parameters under the external magnetic field. We use 
the model Hamiltonian of Shiina et al. within the mean 
field approximation (MFA) for describing the initial state 
of the RXS process in the AFQ phase. 



2 



Jun-ichi IGARASHI and Tatsuya Nagao 




Fig. 1. (a) Antiferro-quadrupolar structure of phase II. Open cir- 
cles are cerium atoms with A and B representing sublattices. 
Solid small circles are boron atoms, (b) Magnetic structure of 
phase III predicted by the neutron scattering experiment. Only 
Ce atoms are shown. The arrows indicate the direction of the 
magnetic moment. The attached numbers represent sites of typ- 
ical directions. 



For T < Tn (= 2.4 K), there appears a magnetic long- 
rangejarder. According to the neutron scattering experi- 
ment ,E3 the ordering pattern of the magnetic moment is 
as shown in Fig. 1 (b) . We construct the wave function in 
the magnetic ground state to be consistent with the or- 
dering pattern, and use it as the initial state of the RXS 
process. Our calculation is limited to T = 0, since an 
intersite interaction correctly reproducing the magnetic 
phase has not been derived yet. 

In the intermediate state of the Ei process, the 5d 
states of Ce are involved so that they have to be mod- 
ulated in accordance with the superlattice spots. Since 
the 4/ states are so localized in space that their coupling 
to lattice is very small. Actually the lattice distortion 
associated with the AFQ order has not been observed. 
Therefore, it is highly possible that the modulation is 
brought about by the Coulomb interaction between the 
5d states and the orbital ordering 4/ states. Introducing 
a reasonable density of states (DOS) for the 5d states, 
we solve a scattering problem of the photo-excited 5d 
electron; the scatterer is a complex of a 2p hole and a 4/ 
electron. 

Using the solution in the intermediate state, and com- 
bining it to the result in the initial state, we calculate the 
RXS spectra on an AFQ superlattice spot G — (5^5) 
and on a magnetic superlattice spot G = {jj^) in the 
magnetic ground state. In the AFQ phase, we obtain 
sufficient intensities on the AFQ spot without assum- 



ing any lattice distortions, thereby demonstrating the 
mechanism that the RXS intensity is brought about by 
the modulation of the 5d states in accordance with the 
AFQ order through the intra-atomic Coulomb interac- 
tion. This situation is different from that of transition- 
metal compounds, where there exists a sizable lattice dis- 
tortion which is the primary origin of the RXS spectra. 
The calculated temperature and magnetic field depen- 
dences reproduce well the experiment .0* We find that 
the azimuthal-angle dependence is closely related to the 
symmetry of the AFQ order, independent of the details 
of model. In the magnetic ground state, we find the RXS 
spectra on the AFQ superlattice spot as a smooth exten- 
sion from the AFQ phase, suggesting that the magnetic 
order has little influence on this spot. We find a small 
pre-edge peak in the electric quadrupolar {E2) process 
thereEj^ On the magnetic superlattice spot, we obtain a 
finite intensity, but it is about two order of magnitude 
smaller than that on the AFQ spot. 

This paper is organized as follows. In Sec. II, we for- 
mulate the RXS spectra as a second order optical pro- 
cess. In Sec. Ill and IV, we describe the initial and the 
intemediate states of RXS. In Sec. V, we discuss the cal- 
culated RXS spectra in comparison with experiments. 
Section VI is devoted to concluding remarks. 

§2. Formulation for Resonant X-Ray Scattering 

The conventional RXS geometry is shown in Fig. 2; 
photon with frequency to, momentum and polarization 
^ (— a 01 tt) is scattered into the state with momentum 
k/ and polarization /i' (= a' or tt'). The scattering vector 
is defined by G = k/ — k^. The cross section Jbr the 
scattering vector G consists of three termsiO'EZlEa' 




Sample 



Fig. 2. Scattering geometry of x-ray scattering. Incident photon 
with wave vector k^ and polarization or tt is scattered into the 
state with wave vector kj and polarization a' or tt' at Bragg 
angle 8. The sample crystal is rotated by azimuthal angle ip 
around the scattering vector G = ky — k^. 
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exp(— iG • Tj) 



-L(G,j).A" + S(G,j)-B 



(2.2) 



Mo,a,{G,uj) = ^^p„(j)exp(-iG • Vj)mu}l^ 



hu! — {Ea 

1 ^fc2 



iT 



En{j)) 



(2.3) 



^ ma;i,„(V'n(j)|z^(j)|A')(A'|zy(j)|V'nU)) 
;ia;-(-BA'-£;„(j)) + ir 

Here c is the velocity of photon and m is the electron 
mass. Note that the cross section is an order of the 
number of Ce sites. 

The first term (eq. (2.2)) represents a non-resonant 
term for the magnetic superlattice spots. Since the mag- 
netic moment comes mainly from 4/ states which are 
well localized, it may be a good approximation to assign 
the moment to each site. The L(G,j) and S(G,j) are 
the form factors of the orbitaLapd spin angular momenta 
at site j, which are given byES'B 

L(G,j) = i5]p„(j) 



X (V'„(j)|/(-G • r)£ + t,f{-G ■ r)|V'„(j)I?.5) 



S(G,j) 
with 



P«(j)(V'«(j)|e-'^-''s|^„(j)), 



fix) = 2 



OO 

y 



(ixY 



2)m\' 



(2.6) 



(2.7) 



Here \tpnij)) represents the initial state at site j, which 
will be evaluated within the MFA discussed in the next 
section. We take the thermal average with probability 
Pn{j)- Operators £ and s represent the orbital and spin 
angular momenta with the center of site j. The L(G, j) 
and S(G, j) converge to the local orbital momentum and 
the spin angular momentum with G — > 0. The scattering 
geometry is contained in quantities A" and B, which are 
defined by 

A" = A' - (A' • G)G, A' = -4sin2 6'(e' x e),(2.8) 



B = e X e+{kf X e){kf ■ e) 
- (k/ X e') X (ki X e), 



(k, X e)(k, 



(2.9) 



where e and e are the initial and scattered polarizations, 
and ki , k/ , and G are normalized vectors of k^ , k/ , and 
G. 

The second term in eq. (2.1) describes a resonant term 
by the Ei process, where an electron in 2p states is virtu- 
ally excited to 5d states and subsequently is recombined 
with the core hole. Since the 2p states are well localized 
around Ce sites, it is a good approximation to describe 



the scattering tensor as a sum of the contribution from 
each site of the core hole. The initial state |V'n(j)) at 
site j has an energy En{j). The intermediate state |A) 
consists of an excited electron on 5d states and a hole 
on 2p states with energy Ej^. = (-E'a — En)/h. The 
life-time broadening width F of the core hole is assumed 
to be 2 eV. The dipole operators XaiJYs are defined as 
xi{j) — X, X2{j) — y, and 2:3 (j) = z in the coordi- 
nate frame fixed to the crystal axes with the origin lo- 
cated at the center of site j. The dipole matrix element 
Adp = {5d\r\2p) = R5dir)rR2p{r)r'^dr is implicitly 
included as a square in the expression (Rsdir) and R2p{r) 
are the radial wave functions for the 5c? and 2p states, 
respectively). It is estimated as Adp = 3.67 x 10~^^ cm 
for a Cej£+ atom within the Hartree-Fork (HF) approx- 
imationO The and F'^ are geometrical factors for 
the incident and scattered photons, respectively, which 
are explicitly written in the Appendix of ref. 31. 

The third term in eq. (2.1) describes a resonant term 
by the E2 process, where an electron in 2p states is virtu- 
ally excited to 4/ states and subsequently is recombined 
with the core hole. In eq. (2.4), k is the wavenumber of 
the incident (and scattered) photon, which is ~ 2.91 x 10* 
cm^^ around the Lm edge. The intermediate states |A') 
consist of an excited electron on 4/ states and a hole 
in 2p states with energy E\i . Quadrupole operators are 
defined as zi = (V3/2)(a;2 _ ^2^)^ ^ {l/2){3z^ - r^), 
Z3 = VSyz, Z4 = ^/3zx, and 25 = VSxy in the coordinate 
frame fixed to the crystal axes. The quadrupole matrix 
element Afp = {Af\r^\2p) = R4fir)r^R2p{r)r^dr is 
included as a square in the expression {R4f{r) is the 
4/ radial wave function). It is estimated as Afp = 
5.69 X l(t,'^° cm^ for a Ce'^+ atom within the HF approx- 
imationEj'' The Q'^ and Q'^^ are geometrical factors for 
the incident and scattered photons, respectively, which 
are explicitly written in the Appendix of ref. 31. 



§3. Initial State 

Each Ce atom 
configuration, '^F; 



5/2 



is approximately in the 4/^- 
in CcBg. The 4/ states are so lo- 
calized that their wave functions are well described to 
be atomic in the HF approximation. The cubic crys- 
tal field lifts the degeneracy; quadruply degenerate Fg 
states have a lower energy than the doubly degenerate 
F7 states. Since the Fg-Fy separation energy is estimated 
as large as ^ 540 K, it is sufficient to consider only the 
Fg states in the description of the initial state, which are 
explicitly written as 



l) = 
T) = 
i) = 



(3.1) 



The state \AI) has the z-component M of the total an- 
gular momentum, where the z axis is taken along the 
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[0, 0, 1] direction. Symbols t{— ±) and cr(=t, i) in \t, cf) 
represent non-Kramers' and Kramers' pairs, respectively. 

3.1 Quadrupolar Ordering Phase (T^ <T < Tq) 

Some quadrupolar ordering is established (Phase II) 
below Tq = 3.2 K. An effective Hamiltonian, which is 
derived on the,basis of the RKKY interaction, is known 
to work well:E3 

H = dJ2[{^ + • + rfrj + err, ■ a, 

+^(^^^;■+^, -^j+c-c,) 

+g^lBY.J^■H, (3.2) 

i 

with 

J=^(t + ^77), (3.3) 

fj. = (2r2'cr^, 2r2'CTy, 2T^cr^), (3.4) 

T'-(r^T-), (3.5) 

J7 - ((V3r^ - T^)c^^ (-%/3t^ - r^)cr^ 2tV^),(3.6) 

C = H-VSt' - T^)c^^ (\/3t^ - r^)cr^ 2tV^),(3.7) 

where represents the sum over nearest neighboring 
Ce pairs. Operators t and a represent the spin matrix 
acting on the variables t and a of the state |tct), respec- 
tively. This system has been p-ytprisivply studied within 
the MFA by Shiina et al;B00B S ~ 0.2 and e - 1 
are known to be suitable for CeBg. In the following, we 
simply summarize the result of the MFA in connection 
to the RXS spectra. 

In the absence of the magnetic field, an AFQ order 
is set in, as shown in Fig. 1(a). We have three types 
of possible ordered phase, in which one of the stag- 
gered quadrupole moments, (Oxy) (= 4(tj,CTz)), (Oyz) 
(= 4:{Tyax)), and (Ozx) (= 4:{Tyay)), is finite. Here {X) 
indicates the thermal average of operator X. We simply 
call them as the Oxy, Oyz, and Ozx phases. For exam- 
ple, in the Oxy phase, applying the MFA to eq. (3.2), we 
obtain the eigenstates of the Hamiltonian at site j: 



\Mj)) = 


71^'^ 


-T)+i|- 


T)}, 


(3.8) 


\Mj)) = 




-i>-i|- 


i>}, 


(3.9) 


\Mj)) - 




-T)-i|- 


T)}, 


(3.10) 


\Mj)) = 


71^'^ 


-i>+i|- 


i>}, 


(3.11) 



with the eigenvalues, 

EiU) = E2{j) = TzD{l + S)\A{Tyaz)l (3-12) 

EsU) = E4j) - ±zDil + S)\A{Tyaz)l (3.13) 

where the upper(lower) sign in eqs. (3.12) and (3.13) is 
for A(B) sublattice. The A{Tyaz) is the staggered or- 
der parameter self-consistently determined, and z (= 6) 




is the number of nearest neighboring pairs. Note that 
Kramers' pairs, and \^2), as well as and 

are still degenerate in Phase II. The probabil- 
ity Pnii) appeared in the preceding section is given by 
cx exp(-i;„(j)/r). Note that pi{j){= P2{j)) > P3(j)(= 
P4{j)) at A sublattice, and vice versa at B sublattice. 

The angular dependences of the charge density for 
those states are defined by 

^''^^^'^^^ 2\R,]ir)\^ E \ir,0,^\4'n)\^ (3.14) 

We calculate these quantities Jising the atomic func- 
tion in the HF approximation.!^ As shown in Fig. 3, 
C^'^{6, (j)) is along the [1, -1, 0] direction, while C^ '^iO, (j)) 
is along the [1,1,0] direction. The charge distribution af- 
ter thermal average is along the [1, —1, 0] direction at A 
sublattice, since Pi(i)(= P2U)) > P3{j){= Pi{i)) there. 
Such anisotropy leads to a modulation in the Sd states 
through the Coulomb interaction in the RXS process. 

The magnetic field induces the staggered octupole mo- 
ment in addition to the staggered quadrupole moments, 
as shown in the upper panel of Fig. 10. As shown in the 
same figure, the octupole moment has little influence on 
the RXS spectra. The magnetic field also induces the 
uniform dipole moment. Recently Saitoh et al. have car- 
ried out a neutron diffraction experiment under magnetic 
field, and have reported the magnetic form factors.E3 
They have argued that a considerable amount of the 
magnetic moment is distributed around B atoms. In this 
context, it may be instructive to evaluate the magnetic 
form factors on Ce atoms, although they are not directly 
related to the present RXS study. We use the atomic 
wave functiDn within the HF approximation in eqs. (2.5) 
and (2.6).t2l' Figure 4(a) shows the form factors on Ce 
atoms for various values of G. iJ = 8 T (H || [0,0, 1]), 
T — 1.5 K. The orbital moment is much larger than the 
spin moment, and the form factors decrease monotoni- 
cally with increasing valupi.|aE |G|. These values do not 
fit the experimental data.BEj'' 



Fig. 3. Charge distributions Ci'2{e,0) and C^'4(g 0-)^ 
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(a) Phase II (T=1 .5 K, H=8 T, H//[0,0,1]) 
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(b) Magnetic Ground State 
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Fig. 4. Absolute values of the orbital moment form factor, 
\L{G,j)\, and of the spin moment form factor, |S{G,j)|, as a 
function of |G|{= 47rsin6»/A). (a) Phase II (T = 1.5 K, H = 8 
T, H II [0,0,1]. Thick solid and thick broken lines represent 
\L(G,j)\ for site j belonging to A and B sublattices, respec- 
tively. Thin solid and thin broken lines are for |S(G,j)|. (b) 
Magnetic ground state. Thick solid and thick broken lines rep- 
resent \L{G,j) \ for site j belonging to site 1 and 3, and belonging 
to site 2 and 4, respectively. Thin solid and thin broken lines are 
for \S{G,j)\. 



3.2 Magnetic Phase (T <T^) 

With further decreasing temperature (T < Tn), a 
magnetic long-range order appears (Phase III). The or- 
dering pattern predicted by the neutron scattering ex- 
periment is shown in Fig. 1(b). The magnetic moment 
m(rj ) is directing on the ab plane; it is given by 



m(rj) cx I cos(^ki • + ^ 



+ 



(' 



cos I k'l 



cos 



TT 

r," 

J 4 



2 • r,- 



"kj' 
(3.16) 



where uj^^ and uj^.^ are unit vectors along the [—1, 1,0] 
direction and along the [1,1,0] direction, respectively, 
andki = {\\\), k2 = (Hi), k'l = (iiO), k'2 = (i|0). 
Of course, there must exist other magnetic domains in 
which the moments are directing on the he and ea planes. 

The magnetic domain described by eq. (3.16) is ex- 
pected to come from the Oxy phase with splitting the 
degeneracy of Kramers' doublet. We have no reliable ef- 
fective inter-site interaction between Kramers' doublets; 
the inter-site interaction given by eq. (3.2) cannot de- 
scribe the magnetic state. Therefore we are satisfied to 
derive the ground state wave function consistent with the 
distribution of the magnetic moment experimentally de- 
termined. First we note that the angular momentum op- 
erator J along n (= (cos (p, sin 0, 0)) is represented within 



the space of jV'i) and |'02) as 

n.J = A|V^i)(V^2|+A*|V^2)(V'i|, (3.17) 

with A — — (1/3) exp(— 10) — (i/\/3) exp(i(/)), while it is 
represented within the space of jV'a) and \^J4) as 

n • J = S|V^3) {M + B* 1^-4) (^-31, (3.18) 

with B = —(1/3) exp(— i0) + (i/V3) exp(i0). Considering 
the magnetic ordering pattern in Fig. 1(b), and noting 
that J is pointing to the direction opposite to the local 
magnetic moment vector, we seek the eigenstate with the 
negati ve eig envalue for eq. (|3.17| ) with (j) = 37 r/4 fo r site 
1, eq. (|3.18i) with </> = -37r/ 4 for site 2, eq. (|3.17| ) with 
(j) = — 7r/4 for site 3, and eq. ( ^.18 ) with = 7r/4 for site 
4. Thus we have the wave function in the ground state 
as 

-i37r/4| 




HV'2)), 
1^4)) , 

+ IV'4)) , 



for site 1, 
for site 2, 
for site 3, 
for site 4, 



(3.19) 

with the eigenvalue —0.91. Considering the g-factor 6/7, 
we have the local magnetic moment 0.78/iB, which is 
close to the value Ojfifi/iB from the analysis of the ^^B- 
NMR measurement .lS 

The orbital and spin form factors are evaluated by us- 
ing the above wave functions. Figure 4(b) shows the 
calculated result as a function of |G|. The form fac- 
tors decreases monotonically with increasing values of 
|G|, which behavior is similar to that in Phase II under 
the magnetic field. We need these quantities to evaluate 
the non-resonant term, eq. (2.2), in the RXS spectra. 
However, as shown later, the non-resonant term gives 
much smaller contribution than the resonant terms in 
the ground state. 

§4. Intermediate State 

In the El process, an electron is excited from 2p states 
to 5d states at a Ce site. The 2p-core hole states are 
split into the states with jp = 3/2 and jp = 1/2 {jp is 
the total angular momentum) due to the strong spin- 
orbit interaction. In the following, we consider only the 
jp — 3/2 states (Lm edge). 

Different from the 4/ states, the 5d states are rather 
extended in space, so that they form an energy band with 
width ~ 15 eV through the hybridization with boron p 
states. We use a shape showii.in Fig. 5 as a model of the 
5c? density of states (DOS). £3 We disregard the depen- 
dence on the symmetries xy, yz, zx, x^ — y^, and 3z^ — r^. 
It is assumed to be occupied by one electron per Ce site 
in the initial state. Such a rather arbitrary choice of the 
5c?-D0S may be justified in a semi-quantitative study, 
since the RXS spectra is not sensitive to the shape and 
the filling of the 5d-D0S. The retarded Green's function 
for the excited 5d electron is defined by 



huj 



+ i5 



rde. 



(4.1) 



where p (e) is the 5(i-D0S, and ep is the Fermi en- 
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Center Center Neighbor 

Fig. 5. Schematic view of the configuration of the intermediate 
state of the Ei process: (a) the 5d excited electron is at the 
central site where the core hole is sitting; (b) the 5d excited 
electron is at the neighboring sites. The 5(i-DOS is schematically 
shown on the upper part. Electrons are occupied below the Fermi 
level. The wavy lines represent the Coulomb interaction taken 
into account. Arrows indicate spins of 5d states. 



Table I. Slater integrals and the spin-orbit interaction for Ce' 
atoms in the Hartree-Fock approximation (in units of eV). 



F'=(4/,4/) 


F'=(2p,5d) 


F'=(2p,4/) 


F'=(4/, 5d) 


F° 26.08 


15.58 


F° 37.68 


F° 13.92 


12.43 


F2 0.568 


1.540 


F^ 3.810 


7.807 






F* 1.894 


5.618 










G'=(2p,5d) 


G'=(2p,4/) 


G'^(4/, 5d) 




Gi 0.485 


G2 0.144 


Gi 1.633 




0.286 


G-* 0.093 


G^ 1.397 








G5 1.086 


(4^ = 0.132 


Csd = 0.138 









*In the RXS calculation, the above values of the anisotropic terms 
are reduced by multiplying a factor 0.8, while the values for 
F^{nl, n'V) are replaced by much smaller values, F^{Af, 5d) = 3.0, 
F0{4/,4/) = 7.0, F''(2p,5d) = 4.0 and F°{2p,4f) = 12.0. 



Note that the energy of 5d states included im- 
phcitly an average interaction with electrons in 4/ states. 

Now we consider the resolvent l/(?iti> — i?int) with Hint 
being the Hamiltonian in the configuration that there are 
one excited electron in the 5d band, one electron per Ce 
site in the 4/ states, and one hole in 2p states at the cen- 
tral site (Fig. 5(a)). We neglect the screening effect by 
the already occupying electrons in the 5d band. When 
the excited "5d electron" is away from the central site, it 



can freely move; the Green's function, eq. (4.1), describes 
this motion. At the same time, the "4/ electron" is inter- 
acting with the "2p hole" on the central site (Fig. 5(b)). 
First, considering the Coulomb interaction between the 
4/ electron and the 2p hole and also the spin-orbit inter- 
action at the central site, we solve the eigenvalue problem 
for the complex of the 4/ electron and the 2p hole. In 
this calculation, we use the Slater integrals evaluated by 
the HF approximation in a Ce'^"'' atom, which are listed 
in Table l& Let |A) be the eigenstate with energy Ex. 
On the other hand, when the 5d electron comes onto the 
central site, it interacts with the 2p hole as well as the 4/ 
electron. We treat this system as a scattering problem 
of the 5d electron, in which the scatterer has 56 (4 x 14) 
internal degrees of freedom specified by A on the central 
site. Taking account of the multiple scattering of the 5d 
electron by this scatterer, we obtain the expression of the 
resolvent at the central site as 

1 



hui — Hint + iS , 
— [G'^'^(?iW + ir — Ex) ^ 5xX' ^ra'^-m"^^ s'^-s"^ 
~ ^m''- s''- \;m"^ s"^ \'\ : 



(4.2) 



where [mf's'^) specifies a state of the 5(i electron. Ma- 
trix Vjyidgdx-.m'ds'dx' includes the Coulomb interaction be- 
tween the 5(i electron and the 2p hole and between the bd 
electron and the 4/ electron. We have to exclude the av- 
erage interaction energy between the 5(i and 4/ electrons 
from this matrix, since it is already included into the en- 



ergy of the bd band. The right hand side of eq. (4.2) is a 
matrix with dimensions 560 x 560 (560 = 10 x 56), which 
we numerically invert. 

Once we obtain the resolvent, we can calculate the 
scattering amphtude, eq. (2.3), by using the relation, 

|A)(A| 



^ huo-{Ef,-En{3))+iT 



E E 



huj Hint + i<5/ m'is'i\;m"is"i\' 



(m"'s"^A'| 



(4.3) 



for j at the central site. It should be noted here that this 
resolvent is the same at all sites of core hole. The scatter- 
ing amplitudes become different at different sublattices 
after multiplying the matrix elements of the dipole op- 
erators between the initial and the intermediate states. 
Using the resolvent, we can also calculate the absorption 
coefficient A{lo) in the Ei process. 



A{uo) « ^^p„(j)(^„(j)|2:"|mVA) 



J na 



X ( — — I Im 



1 



■kJ \^ - Hint+iSj ^d,dx-m'ds'dX' 

X {m"'s"'\'\x"\Mj)), (4.4) 

where ImX indicates the imaginary part of the quantity 
X. 

In the E2 process, an electron is excited from 2p states 
to 4/ states at Ce sites. Two electrons occupy the 4/ 
states on the Ce site. Since the 4/ states are well local- 
ized, it may be sufficient to consider only the core-hole 
site for the intermediate state |A') in eq. (2.4). We ob- 
tain I A') by numerically diagonalizing the Hamiltonian 
matrix within the space of two 4/ electrons and one 2p 
hole. We fully take account of the multiplets, where the 
necessary Slater integrals as well as the-|Spin-orbit inter- 
action parameter are listed in Table lS3> 
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Fig. 6. Absorption coefficient A{ui) of the Ce Lni edge in com- 
parison with the experimental spectra (ref. 37). T = 2.7 K and 
H = 0. Solid and broken lines represent the spectra calculated 
with F^{2p,5d) = 3 eV and 0, respectively. 



§5. RXS Spectra 

5.1 Quadrupolar Ordering Phase (T^ <T < Tq) 

Before going to the discussion of the RXS spectra, we 
first calculate the absorption coefficient A{ui) with the 
help of eq. (4.4). The contribution from the E2 process 
is safely neglected, since it is too small by the estimate 
of the HF transition matrix elements. Figure 6 shows 
the calculated result at T ^2.7 K, H — Q, in com- 
parison with the experiment .EZP The core- hole energy is 
adjusted such that the peak is located at huj = 5722 
eV. Note that the temperature dependence is negligible. 
The spectrum is a reflection of the Sd-DOS; the attrac- 
tive interaction between the 5d electron and the core hole 
makes the peak move to the low-energy region. The cal- 
culated spectrum corresponds well to the experimental 
shape around the Lm edge, indicating the appropriate- 
ness of the assumed 5d DOS. Considerable intensities in 
the high energy region in the experimental spectra may 
come from the mixing of 5c? states to other states such 
as 6s states of Ce and 3s states of B.Now we calculate 
the RXS spectra for G = (|^^), following the procedure 
in the preceding section. Figure 7 shows the calculated 
spectra as a function—pf photon energy, in comparison 
with the experiment .EJ r = 2.7Kandi7 = 0. The az- 
imuthal angle ip is set to be zero such that the scattering 
plane contains the [1,-1,0] crystal axis. The relative 
volumes of three domains, the Oxy^ Oyz and Ozx phases, 
are not known in the experiment, so that we have tenta- 
tively averaged the contributions from three phases with 
equal weight. Since the dependence on the photon en- 
ergy is the same in three domains, the spectral shape 
is not influenced by the change of the relative volume 
of domains. Only changeable are the relative intensities 



between the a a' channel and the cr — > tt' channel; 



in the O 



xy 



phase, the intensity of the a 



tt' channel 

is larger than that of the a — > cr' channel. Such polar- 
ization analysis has not been done in the experiment. L3-' 
We obtain an one-peak structure from th£ E\ process, 
in good agreement with the experiment .0 The contri- 
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Fig. 7. RXS spectra for G = (^^^) in comparison with the ex- 
periment (ref. 14). T = 2.7 K, H = 0, and ip = 0. Intensities 
from three domains are averaged with equal weight. The solid 
and broken lines represent the intensities for the a —> a' channel 
and the a ^ n' channel, respectively. 




Fig. 



RXS intensity of the peak at hw = 5722 eV for G 







(^ii) in the Oxy phase, as a function of temperature. H 
and i/j = 0. The solid and broken lines represent the intensities 
for the a ^ a' channel and the it — » 7r' channel, respectively. 
Solid circles at T = is the peak intensity evaluated in the 
magnetic ground state. The upper panel shows the staggered 
quadrupole moment (Oxy) and its square {Oxyf in the Oxy 
phase within the MFA. 



bution of the E2 process is two order of magnitude less 
than that of the Ei process, so that we have no pre-edge 
peak visible. 

As mentioned before, the intensity arises from the 5d 
states modulated by the anisotropy in the 4/ charge dis- 
tribution through the intra-atomic 5d-4/ Coulomb in- 
teraction in the intermediate state. Figure 8 shows the 
intensity of the peak at hw = 5722 eV in the O^y phase, 
as a function of temperature. The curves are extended 
to r < Tn by assuming the Oxy phase. Its temperature 
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Fig. 9. RXS intensity of the peak at hu) = 5722 eV for G = 
(iii), as a function of azimuthal angle ■(/>. T = 2.7 K and 
H = 0. (a) the tr — ► cr' channel; (b) the a ^ n' channel. The 
solid, broken, and dotted lines represent the contributions from 



the Ox 



phases, respectively. 



dependence seems similar to that of (Oxy)'^, indicating a 
direct reflection of the AFQ order. 

Another important quantity is the dependence on the 
azimuthal angle. Figure 9 shows the RXS intensity of 
the peak at hoj — 5722 eV. The contributions from three 
domains are separately shown. The dependence of the 
a —f a' channel is quite different from that of the ct — > tt' 
channel. The curves in the Oxy, Oyz and Ozx phases 
can be transformed into those for the Oyz, O^xi and Oxy 
phases, respectively, by shifting ip with 27r/3. This three- 
fold symmetry around G = (555) perfectly matches the 
relation between the order parameters of three domains. 
Thus the azimuthal-angle dependence is closely related 
to the geometry of scattering as well as the symmetry of 
the AFQ order, not related to the details of the model. 
Therefore, the examination of the azimuthal angle de- 
pendence may be useful to determine the symmetry of 
the AFQ phase. 

The initial state is rather sensitive to the magnetic 
field. For example, as shown in the upper panel of 
Fig. 10, a sizable staggered octupole moment {Txyz) 
(= 2(r'')) is induced by applying the magnetic field along 
the [0, 0, 1] direction. The lower panel of Fig. 10 shows 
the intensity of the peak at hu) = 5722 eV as a function 
of magnetic field. The intensity increases only gradu- 
ally with increasing iJ, which behavior is close to the 
variation of the staggered quadrupole moment. This in- 
dicates that the induced staggered octupole moment has 
little influence on the RXS spectra. Finally in this sub- 
section, we demonstrate in Fig. 11 that the calculated 
temperature dependence reproduces well the experiment 
for H II [l,l,-2].llJ 

5.2 Magnetic Ground State 

Using the wave functions given in the preceding sec- 
tion for the magnetic ground state, we calculate the RXS 
intensities for the AFQ superlattice spot G = (555) and 
for the magnetic superlattice spot G = (315)- Figure 
12 shows the calculated spectra as a function of photon 
energy for H = 0. 

For G = (^^^), the spectral shape is close to the one 
in the AFQ phase. The intensity of the main peak is a 
smooth extension from the Oxy phase, as shown in Fig. 8. 
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Fig. 10. RXS intensity of the peak at Tiuj = 5722 eV for G = 
(iii), as a function of magnetic field. H || [0,0, 1], T = 2.7 K 
and 1/1 = 0. The solid and broken lines represent the intensities 
for the a ^ a' channel and the ct — > tt' channel, respectively. 
The upper panel shows the staggered quadrupole moment (solid 
line) and the octupole moment (dotted line). 




Fig. 11. RXS intensity of the peak at fiu = 5722 eV as a function 
of temperature, for various values of H (H || [1,1,-2]). G = 
( ^ i i ) . The curve for H = is not the domain average but the 
limit of ^ 0. 



This indicates that the primary origin is the anisotropic 
charge distribution associated with the AFQ order in the 
4/ states, and that the magnetic order, which lifts the 
degeneracy of Kramers' doublet, has little influence on 
this spot. On the other hand, a pre-edge peak around 
fiuj = 5710 eV is enhanced to become visible in the E2 
process. 
For G 



> 4 4 2 ^ 



we have the intensity two order of 
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Fig. 12. RXS spectra in the magnetic ground state as a function 
of photon energy: (a) G = (b) G = (jj^)- = and 

■0 = 0. The pre-edge peak is magnified in the inset. The solid 
and broken lines represent the intensities for the ct — > cr' channel 
and the cr — ► tt' channel, respectively. 

magnitude smaller than for G = (^55)- The origin is 
purely magnetic, that is, the main peak arises from the 
modulation of the 5d states in accordance with the mag- 
netic order through the intra-atomic Coulomb interac- 
tion with the 4/ states. The contributions from mod- 
ulation in accordance with the AFQ order are canceled 
out among different magnetic sublattices. No intensity 
comes out in the cr — > cr' channel. Note that the contri- 
bution of the non-resonant term (eq. (2.2)) is found too 
small to be seen. 

§6. Concluding Remarks 

We have studied the RXS spectra near the Ce Lm ab- 
sorption edge in the AFQ phase and also in the magnetic 
ground state, on the basis of a microscopic model that 
the 2p and 4/ states of Ce sites are atomic while the 
5c? states form an energy band with a reasonable DOS. 
Since the lattice distortion associated with the AFQ or- 
der has not been observed, we have paid no attention on 
it in the present calculation. The Fs states of the Ce 4/ 
states have a lower energy than the F7 states under the 
crystal field. 

In the AFQ phase, we have employed an effective inter- 
site interaction between the Fg states, which is developed 
by Shiina et. ai, and applied the MFA for describing 
the initial state. In the description of the intermedi- 
ate state, on the other hand, those energies can be ne- 
glected in comparison with much larger energies such as 
the intra-atomic Coulomb interaction and the energy of 
the 5c? band. Introducing the Green's function, and tak- 
ing full account of the intra-atomic Coulomb interaction, 
we have solved the scattering problem of the excited 5c? 
electron by the core-hole potential. We have found that 
the initial state has a sensitive effect on the RXS spectra 
through the dipole matrix elements. 



We have obtained relatively large RXS intensities on 
the AFQ superlattice spot without assuming the lattice 
distortion, thereby demonstrating the mechanism of the 
Coulomb interaction. This situation is different from 
that of transition- metal compounds. The temperature 
and magnetic field dependences of the RXS spectra re- 
produce well the experiment Jlil' We have found that the 
azimuthal-angle dependence is closely related to the sym- 
metry of the AFQ order. Therefore this quantity may be 
useful to determine the order parameter. We hope this 
study prompts experimentalists to measure this quantity. 

In the magnetic ground state, we have constructed the 
wave function to be consistent with the ordering pattern 
of the magnetic moment determined by the neutron scat- 
tering experiment. Using this as the initial state, we have 
calculated the RXS spectra on an AFQ superlattice spot 
and on a magnetic superlattice spot. The intensity on the 
AFQ spot is a smooth extension from the AFQ phase. In 
addition, we have found a small pre-edge peak in the E2 
process. Since the main peak intensity will be reduced 
by the absorption correction, the pre-edge peak might be 
observed in future experiment. On the magnetic super- 
lattice spot, we have a finite but much smaller intensity 
than that on the AFQ spot. It may be hard to confirm 
experimentally the spectra on this spot. 

We have assumed the shape of the 5c?-D0S rather ar- 
bitrarily. This is partly justified by the fact that the 
characteristics of the RXS spectra discussed in this pa- 
per do not sensitively depend on the details of the 5d 
DOS. Nevertheless, a band structure calculation for the 
5d states may be necessary for more quantitative study. 
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